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Abstract 

We consider multiple orthogonal polynomials corresponding to two Macdonald 
functions (modified Bessel functions of the second kind), with emphasis on the 
polynomials on the diagonal of the Hermite-Pade table. We give some properties 
of these polynomials: differential properties, a Rodrigues type formula and explicit 
formulas for the third order linear recurrence relation. 



1 Macdonald functions 

In this paper we will investigate certain polynomials satisfying orthogonality properties 
with respect to weight functions related to Macdonald functions K u (z). These Macdonald 
functions are modified Bessel functions of the second kind and satisfy the differential 
equation 

z 2 u" + zv! - (z 2 + v 2 )u = 0, 

for which they are the solution that remains bounded as z tends to infinity on the real 
line. They can be given by the following integral representations 

(see, e.g., The Bateman Manuscript Project || Chapter VII], Nikiforov and Uvarov 
pp. 223-226]). Useful properties are 

-2K' v {z) = K u ^{z) + K v+1 {z), (1.3) 

-—K v {z) = K u ^(z)-K u+1 {z). (1.4) 

z 
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These functions have the asymptotic behavior 



.2z 

and near the origin 



KAz) = [^-) 1/2 e- z [l + 0(l/z)}, z^oc, (1.5) 



z v K v (z) = 2 u - l T(u) + o(l) if 2 -> 
AToW = -log ^ + 0(1), if z^O 



;i.6) 



The moments of these functions are given explictly by 

r x»K v {x) dx = 2^ 1 r( /i + V - + - )T( - ~ ^ + K(u±z/+1)>0. (1.7) 

Jo 2 2 

We will use as weight functions the scaled Macdonald functions 

p„(x) = 2x u/2 K u {2Vx), x>0. (1.8) 



which behave like exp(— 2a/x) as x — > oo. From the properties ( |1.3| ) and (|1.4 ) one easily 
obtains the differential properties 

(aTW = -aT^+V+i, (1.9) 
pUi = -P- (1-10) 

so that p v satisfies the second order linear differential equation [x u+l (x~ v p u )']' = p v . By 



a simple change of variable in (177) we can find the moments of p v which are given by 



Jo 



x» Pu (x)dx = r(p + v + i)r{p + i), (i.ii) 



in particular we see that the nth moment of po is given by (n!) 2 . 

In || A. P. Prudnikov formulated as an open problem the construction of the orthog- 
onal polynomials associated with the weight function p u (and for more general weights, 
known as ultra-exponential weight functions). In the present paper we will show that it is 
more natural to investigate multiple orthogonal polynomials for two Macdonald weights 
p v and p v+ \ since for these multiple orthogonal polynomials one has nice differential prop- 
erties, a Rodrigues formula, and an explicit recurrence relation. 



2 Multiple orthogonal polynomials 

We will investigate two types of multiple orthogonal polynomials for the system of weights 
(Pu,Pu+i) {y > 0) on [0, oo) with an additional factor x a . Let n, m e N and a > — 1, 
then the type 1 polynomials (A® m , B% m ) are such that A^ m is of degree n and B" m is of 
degree m, and they satisfy the orthogonality conditions 

POO 

/ [A^ m (x)p v (x) + B^ m (x)p v+1 {x)]x k+a dx = 0, k = 0,1,2,... ,n + m. (2.1) 
J o 



2 



This gives n + m + 1 linear homogeneous equations for the n + m + 2 unknown coefficients 
of the polynomials A" m and -B" m , so that we can find the type 1 polynomials up to a 
multiplicative factor, which we will fix later. It will be convenient to use the notation 

Qn,m( X ) = A-n, m {x)Pv (x) + B™ m {x)p v+ \ (x) . 

Type 2 polynomials p% m are sucn that p% m is of degree n + m and satisfies the multiple 
orthogonality conditions 

P%, m (x)p 1 ,(x)x k+a dx = k = 0,1,2,... ,n- 1, (2.2) 

Pnm( x )pu+i( x ) xk+a dx = k = 0,1,2,... ,m- 1. (2.3) 





This means that we distribute the n + m orthogonality conditions over the two weights 
x a p u and x a p u+ i. This gives n+m linear homogeneous equations for the n+m+1 unknown 
coefficients of p^m- For type 2 polynomials we will consider monic polynomials, thereby 
fixing the leading coefficient to be 1 and leaving n + m coefficients to be determined from 

(PD-(PD. 

Multiple orthogonal polynomials are related to Hermite-Pade simultaneous rational 
approximation of a system of Markov functions near infinity, (see, e.g., Nikishin and 
Sorokin |j] and Aptekarev [0]). In the present situation the functions to be approximated 
are 

.In Z X In Z X 



Type 1 Hermite-Pade approximation (Latin type) consists of finding polynomials A 



n,mj 



B Ujm , and C n ^ m such that 



Avn(*)/lCz) + B n;m ( Z )f 2 (z) - C rhm (z) = 0(z n m ') Z^OO, 

and type 2 Hermite-Pade approximation (German type) is simultaneous rational approx- 
imation of (fi, f 2 ) with a common denominator p n , m , i-e., one wants polynomials p n ,m, 
Rn, m and S njm such that 

Pn,m(z)fl(z) - R n , m ( z ) = 0{z~ n ~ 1 ), Z -> OO, 



Pn,m(z)f 2 {z) - S n . m (z) = 0(Z ™ ), 



z — > OO. 



The polynomials v4 n m , 5„ jm then are precisely the type 1 polynomials m and -B" m and 
the polynomial p n , m is the type 2 polynomial m . The numerator C„ im is then given by 



Cn,m(^) 







M^) h fr+iW 



Z — X Z — X 

and for type 2 approximation the numerators are 

p ( \ _ f°° Pn,m( Z ) ~ Pn,m( X ) / \ a , 
^Mi.m\Z) — / Pu\ x ) x <3X, 

Vn Z - Z 



x a <ix, 



Pn,m\ Z ) Pn. 



X 



£ — X 



Observe that we can write the system (fx, f 2 ) as 



p / \ f°° ti \Pv+i{x)x a f°° p u+ i(x)x a 
fi( z ) = f( x ) dx, f 2 (z) = / dx, 

Jo Z X Jq z X 

where / is itself a Markov function 

= Pv {x) i_ r° s^ds 

J[X) Ph-i(z) Wo (x + s)[JZ +1 (2yrs) + yJVi(W 

which follows from a result by Ismail |§. Note however that / is a Markov function of 
a positive measure for which not all the moments exist. Nevertheless we can still say 
that (fx, f 2 ) is a Nikishin system, which guarantees that the polynomials A" m , B" m and 
Pn,m a U can be computed and their degrees are exactly n, m, and n + m respectively. 
Furthermore the zeros of q™ n , q% + i n , p",n an d Pn+i,n wm an be on (0, oo). 



3 Differential properties 

In this section we will give some differential properties for the type 1 and type 2 multiple 
orthogonal polynomials. We will only consider the multiple orthogonal polynomials on 
the diagonal or close to the diagonal, i.e., when n = m or n = m + l. These are the most 
natural and they satisfy interesting relations. 

Theorem 1 For the type 2 multiple orthogonal polynomials we have for every a > — 1 



^p£„(s) = ZnpSS-iC*), = (2n - Ijp^-iW- (3.1) 



This means that the differential operator acting on type 2 multiple orthogonal polynomials 
lowers the degree by one and raises the a by one. This should be compared with the 
corresponding differential property 

^-L«(x) = -L^l(x) 
ax 

for Laguerre polynomials (see, e.g., Szego |7], Eq. (5.1.14) on p. 102]). The normalizing 
constant is different here since Laguerre polynomials are not monic polynomials. 
Proof: We begin by using (ET 



p^ n {x)p v (x)x k+a dx = 0, k — 0,1,... — 
and then use Q1.10Q to find 

POO 

- / p° n (x)p' v+1 (x)x k+a dx = 0, k — 0,1,... ,n — l. 
Jo 

Integration by parts then gives 

POO 

-Pln(*)pnu + l{x)x k+ X + / bn,n(^ + 1 ' Pu+1 (*) *X = 0, k = 0, 1, . . . , 71 - 1. 

JO 
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The integrated terms vanish for k > 1 and a > — 1 because of the behavior of p v+ \ near 
and oo, given in ( |1.5| ) and ( |1.6| ). Now work out the differentiation of the product, then 
we have 



oo 







\pl n (x)]'x k+a p„ +1 (x)dx 

+ (k + a) I p° n {x)x k+a - 1 p v+1 {x)dx = 0, k = 1,2,... ,n- 1. 



oo 



The second integral is zero for k = 1,2,... , n because of the orthogonality condition 
(|2.3|), hence we conclude that 



oo 



[Pn,n(x)]'x k+a Pu+1 (x) dx = } k = 1 , 2, . . . , 71 - 1 , 

^0 

or equivalent ly 

/>oo 

/ [K.„(x)]'x £+a+1 p,+i(x) = 0, £ = 0, 1, . . . , n - 2. (3.2) 
Next we do a similar analysis with ( |2.3| ), in which we use ( |1.9| ) to find 

> 

Pl n (x)x k+a+u+1 [x- v p v (x)\' dx = 0, fc = 0, 1, . . . , n - 1. 
Integration by parts gives 

-K>)* fc+Q+i p^c 

POO 

+ / [p^ n {x)x k+a+u+l ] x~ u Pu {x) dx = 0, fc = 0,1,... ,71-1. 

The integrated terms vanish for k > and a > — 1. Working out the derivative of the 
product then gives 



oo 

lk+a+1 



\plni X )\ X Pvi?) dx 

poo 

+ (k + a + u+1) / pl n (x)x k+a p u (x)dx = 0, k = 0,1,..., ri-1. 

The second integral is zero for k = 0, 1, . . . , n — 1 because of the orthogonality (|2~2|), hence 
we have 



K„(^)]'^ +a+1 P^(^) dx = 0, fc = 0, 1, . . . , n - 1. (3.3) 

'o 

Now [p° n ]' is a polynomial of degree 2n — 1 with leading coefficient 2ti (since we normal- 
ized the type 2 multiple orthogonal polynomials by taking monic polynomials), and by 
( |3.2| ) and ( |3.3| ) it satisfies the orthogonality conditions (|2T2|)-(|2T3|) for the type 2 multiple 
orthogonal polynomial PnV-i- By unicity we therefore have [p" n (x)]' = 2np ( ^_ 1 (x). A 
similar reasoning also gives the result \Pnn~i( x )]' = ~ ^)Pn-l n~i( x ) ■ Note however 
that the analysis does not work when m ^ {n, n — 1}. ■ 

There is a similar differential property for type 1 multiple orthogonal polynomials 
which complements the differential property of the type 2 multiple orthogonal polynomials 
given in the previous theorem. 
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Theorem 2 For the type 1 multiple orthogonal polynomials we have for every a > 

This means that the differential operator acting on x a times the type 1 polynomials 
raises the degree by one and lowers the a by one. This should be compared with the 
corresponding differential property 

^-[e- x x a L a n {x)\ = (n + l)e~ x x^ 1 L a n l\{x) 

for Laguerre polynomials. What we will really prove is that the derivative of x a q® n (x) 
is proportional to x a ~ 1 q"+l n (x) (and similarly for the derivative of x a q% n _ 1 (x) which is 
proportional to The choice of the proportionality factor one in (3.4) will 



fix the normalization which was left unspecified by the homogeneous system of equations 
&■ 

Proof: We begin by using (|2.1| ), which we can write as 

POO 

/ x a <£ (x h+1 ydx = 0, k = 0,1,... ,2n. 
Jo 

Integration by parts then gives 

POO 

<n(*V +a+1 |r- / [x a q a n Jx)]'x k+1 dx = 0, fc = 0,l,...,2n. 

JO 

The integrated terms will vanish for every a > — 1 and k > 0. For the integrand of the 
integral we have 

bXnO*)]' = KK^Puix)] +[x a Bl n {x) Pu+l {x)] 

= ax a - l [Al nPv {x) + B« n (x)p u+1 (x)} 
+ x a [(Al n (x)y Pu (x)+Al n (x)p' u (x) 

+ i B n,n( X )y 'P«+l( X ) + B n,n{x)p' v+X {x)]. 



Now use (PHPDD to find 



[x a ql n {x)]' = x a -\{{a + v)J% t Sx)+<K,Sx))' ~ xB« n (x))p u (x) 
+ (aB« n (x) - Al n {x) + x(B:Jx))') Pi/+1 (x)]. 

Observe that (a + v)A% n (x) + x(A% n (x))' — xB" n (x) is a polynomial of degree at most 
n+1 and aB^ n (x) — A% n (x) + x(B£ n (x))' is a polynomial of degree at most n, so that 
we can write 

[x a €,rXx)}' = x Q - 1 [P n+ i(x)p I ,(x) + Q n {x)p v+1 {x)\ (3.5) 

with polynomials P n+ i and Q n of degree at most n + 1 and n respectively. We therefore 
have 

poo 

x a - 1 [P n+l (x)p u (x) + Q n (x)p u+l (x)}x k+l dx = U, fc = 0,l,... ,2n. 



In addition to this, we also have 

[ar<£ n (x)]'dx=ar<£ tn {x)\~ = 0, 
so that 



POO 

/ lx a Qn,nM'x £ dx = 0, £ = 0,1,... ,2n + l. 
Jo 



In view of (12.11) and (13.51) this means that 



[x a qn,n(x)}' = constant x a 



X). 



Since q" n and g"+| n are only determined up to a normalizing factor, we can choose the 
constant equal to one, thereby fixing the normalization. 

We can repeat the analysis for q% n _i with minor changes. In this case we have 

K<n-i(*)]' = ^[((a + *K»-iO»0 + - x J B„V 1 (x))p ! ,(x) 

+ (aB^x) - A^ix) + x( J B« n _ 1 (x))')p, + i(x)]. 

and since (a + v)A™ n _ x (x) -\-x{A < ^ n _ 1 {x))' — xB^ n _ x {x) is a polynomial of degree at most 
n and aS" n _ 1 (a;) — A% n _ x {x) + x(i?"„_ 1 (x))' a polynomial of degree at most n also, the 
orthogonality 

roo 

[x a q% n _ x (x)]'x i dx = 0, £ = 0,1,... ,2n 



o 

gives the required result. 

Observe that the reasoning does not work for [x a q" m (x)]' when m ^ {n, n — 1}. ■ 

The proof also shows that 

K+U x ) = (u + v)Al n (x)+x[Al n (x)]'-xB« n (x), 

K+l,n(x) = aB^ n (x)-A^ n (x)+x[B^ n (x)]', 

and similarly 

Kni*) = ^ + ^Al n _ 1 (x) + x[Al n _ 1 (x)]'-xBl n _ 1 (x), 

B:-\x) = aB^ n _ x {x) - A a n>n _ x {x) + x[B^ n _ x {x)]' . 

4 Rodrigues formula 

As a consequence of Theorem ^| we have the following Rodrigues formula for the type 1 
multiple orthogonal polynomials. 

Theorem 3 The type 1 multiple orthogonal polynomials can be obtained from 

J2n (1211+1 

' (x 2n+a p„(x)) = *Xn-i(*), — — (x 2n+1+a p„(x)) = x a ql n {x). (4.1) 



dx 2n V- rvVJJ ~ *n,n-lW, 



7 



Proof: By combining the two formulas in ( |3.4j) we get 

d 2 

(„ol a ( \\ _ a-2 a-2 

^2 V i ln,n—l\ u 'J/ x Hn+l,n- 



Iterate this k times to get 

d 2k 



dx 2k 

Choose n = and a = 2k + (3 to find 

d 2k 



{ X Qn,n-l( X )) — X Qn+k,n+k-V 



dx 2h 



Now qo-i(x) = Aq -.\p v {x), where Aq x is a constant, which we will take equal to one to 
fix normalization. Thus, if we replace (3 by a and k by n, then we find the first formula 
in dpi). 

Similarly, we can use 

and differentiate it 2k times, to find 



2k+l 



d lk 



X). 



(r a n a — ^-2k-l a-2k-l 

d x 2k+l \ L ln,n-l\ u ' ) ) x H n +k,n+k 

Choose n = and a = 2k + 1 + (3 to find 

Mk+l 

Replacing (3 by a and by n then gives the required formula. ■ 

The Rodrigues formula allows us to compute the type 1 multiple orthogonal polyno- 
mials explicitly. There is a relationship between type 1 and type 2 polynomials (which is 
not typical for the Macdonald weights but holds in general) and this allows us to compute 
the type 2 multiple polynomials as well. Indeed, we can write 

1n,n \ _ ( A n n B n n \ ( Pu 
Qn,n-lJ \ ">n-l ^n,n-lj \Pv+l 



so that we have 



Aa r>a J \ n OL J \ n 

^n.n-l SD n,'n—\J \Hn,n~lJ \Pv+l, 

Writing the inverse of a matrix as the adjoint matrix divided by the determinant gives 

B n ^n—l B n n \ I (?n,n | \ A a R a — A a R a 1 I ^ v I (A 0\ 

Aa Aa I I „a I ~ [ f± n,n- D n,n-l /1 n,n-l- D n,nJ I _ I • Y 1 -^) 

A n,n-1 A n,n / \<in,n-lj \Pv+\ J 



s 



The polynomial A^ n B^ n _ 1 — A^ n _ 1 B^ n is of degree at most 2n and satisfies 

> 

[A^B^x) - Al n ^{x)Bl n ^)]Pu{x)x k+a dx 



o 



[Bln- X {x)ql n {x) - Bl n {x)ql n ^{x)]x k+a dx = 0, k = 0, 1, . . . , n - 1 
where we have used (|4.2| ) and (|2.1|) . Furthermore we also have 



n 



o 



[Al^Bl^x) - < n _ 1 (x) J B: >n (x)]p, +1 (x)x fc+Q ^ 

-K,n-i{x)ql n {x) + Al n (x)ql n ^(x)]x k+a dx = 0, k = 0, 1, . . . , n - 1. 

Hence we can conclude that 

^„(*)*£*-i(*) - <„_i(o;)^ B (x) = constant K», 
where the constant is the leading coefficient of the polynomial. A similar reasoning using 

(Qn+l,n\ _ f A n+l,n ^+l,n| [ P" \ 
Qn,n ) V A n,n B n,n ) \Pv+l ) ' 

gives 

A n+l,n( X ) B nA X ) ~ A n,n{ X ) B n+l,n{x) = Constant K+l,n(^)- 

5 Recurrence relation 

To simplify the notation we put 

It is known that the sequence P n {x), n = 0,1,2,... satisfies a third order recurrence 
relation of the form 

xP n (x) = P n +i(x) + b n P n (x) + c nJ R„_i(x) + d n P n _ 2 (x), (5.1) 

(see, e.g., || or ||). Explicit formulas for the recurrence coefficients are given in the 
following theorem. 



Theorem 4 The recurrence coefficients in ( \5. 1\ ) are given by 



b n = ( n + o; + l)(3n + a + 2i/)-(a + l)(i/-l) 
c n = n{n + a){n + a + z/)(3n + 2a + v) 

d n = n(n — l)(n + a — l)(n + a)(n + a + v — l)(n + a + v). 
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Proof: We begin with the recurrence coefficients of even index, which are used in the 
recurrence relation 

X Vl, n ( X ) = Pn+l,n( X ) + b 2nPn,n( X ) + OlnP^n-M + ■ (5.2) 

We will write the polynomials m explicitly as 

n+m 



n+m—k 
Jj , 



k=0 

and since we are dealing with monic polynomials, we have 

<m(°) = 1- 

Comparing the coefficient of x 2n in ( |5.2| ) gives 

hn = < n (l) -CLnl 1 ). (5-3) 

hence we need to know the coefficients a" m (l). Comparing the coefficient of x 2n ~ 2 and 
x 2n-3 respectively in (|3.1| ) gives the recurrence 



(2n - l)< n (l) = 2n<£_ 1 (l), (2n - 2)a« n _ 1 (l) = (2n - 1X^(1). 
Combining these relations gives 

^ 4-2 
a n,n(-0 = 7 a n-l,n-l(l)) 

n — 1 

which leads to 

< n (l)=n<+ 2 "- 2 (l). 

In order to obtain an explicit formula, we compute explicitly by solving the system 
of equations 

oo 

(x 2 + a^ 1 {l)x + a1 l {2))x a p v {x)dx = 



o 



oo 



(x 2 + a^ 1 (l)x + a^ 1 (2))x a p u+1 (x)dx = 0. 



o 



Using the moments flL.ll ) this gives 



a? tl (l) = -2(2 + a)(2 + a + i/), (5.4) 

a? fl (2) = (l + a)(2 + a)(l + a + i/)(2 + a + i/). (5.5) 

This gives 

a^ n {l) = -2n{a + 2n){a + 2n + u), (5.6) 
and since 2na° +l n (l) = (2n + l)a"^ 1 (l) this also gives 

< + in(!) = -(2n + l)(a + 2n + l)(a + 2n + i/ + 1). (5.7) 
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Inserting this in (|5.3| ) gives the requested formula for the recurrence coefficient & 2n . 
Next we compare the coefficient of x 2n ~ l in (|5.2|) to find 



C 2 n = < n (2) - <+i,„(2) - &2n< n (l)- 



(5i 



This means that we also need to know a" n (2) and a" +ln (2). Compare the coefficient of 



x 



2n-3 



and x 2n 4 respectively in (|3 . 1|) , then 



(2n - 2)< n (2) = 2n<£_ 1 (2), (2n - 3)a^ n _ 1 (2) = (2n - l)a^i, n _ 1 (2) 



which combined gives 



a n,n(2) 



(2n)(2n-l) ^ a+2 



(2n-2)(2n-3) 



a n-l,n-l(2) 



(2n)(2n- 1) 



a 



a+2n-2 
1,1 



(2). 



Using ( |5.5| ) this gives 

< n (2) = n(2n - l)(a + 2n - l)(a + 2n)(a + 2n + i/ - l)(a + 2n + v) 
and since (2n — n (2) = (2n + l)a"+ 1 (2) this also gives 



(5.9) 



(2) = n(2n + !)(« + 2n)(a + 2n + l)(a + 2n + v){a + 2n + v + 1). (5.10) 



Inserting (|5.9|) , ( |5.10| ), (|5.6|) and the expression for b 2n into (|5.8| ) gives, after some straight- 
forward calculus (or by using Maple) the requested expression for c 2 „. 
Finally, compare the coefficient of x 2n ~ 2 in ( |5.2j) , then 

(5.11) 



d 2n = < n ( 3 ) - <+l,n( 3 ) ~ fo 2n<„(2) - C 2n a" ^1), 



so that we need a" n (3) and a" +ln (3). To this end we compare the coefficient of x 
and x 2n ~~ 5 respectively in (3.1) to find 

(2n - 3)< n (3) = 2n< n _ 1 (3), (2n - 4)a^„ 1 (3) = (2n - lK^^), 

which combined gives 



2n-4 



°n+l,n(3) 



2n(2n + 1) 



<r-r(3) 



(2n + l)(2n)(2n-l)^ a+2w _ 2 



*2,1 



(3). 



(2n-2)(2n-3) ^ 6 

In order to find a 2 1 (3) we will compute p 2 1 explicitly by solving the system of equations 

/"OO 

/ (x 3 + a 2 l 1 (l)x 2 + a% x {2)x + a^))^^) cfo = 

(x 3 + a% x (l)x 2 + a% x (2)x + a 2 l 1 (3))x a+1 p^(a;) rfx 

(x 3 + a^ljx 2 + a 2il (2)x + ^{Z))x a p v+1 [x) dx 

which by using (|1.11|) and some calculus gives 

a^(3) = -(3 + a + i/)(2 + a + + a + i/)(3 + a) (2 + a)(l + a). (5.12) 








11 



Using this gives 
„« _ (2n + l)(2n)(2n-l) 

a ra+l i nl'Jj — g 

(a + 2n + i/ + l)(a + 2n + z/)(a + 2n + v - l)(a + 2n + l)(a + 2n)(a + 2n - 1), (5.13) 
and since (2n — 3)a° n (3) = 2na c ^_ 1 (3) we also have 

_ 2n(2n-l)(2w-2) 

(a + 2n + i/) (a + 2n + v - l)(a + 2n + 1/ - 2)(a + 2n)(a + 2n - l)(a + 2n - 2). (5.14) 



Using ( p.!4| ), ( p.!3| ), (|5.9|), ( [5.7] ) and the formulas for 6 2 n and C2 n in ( |5.11| ) then gives the 



requested expression for d 2 n- 

In a similar way we proceed with the odd indices which appear in the recurrence 
relation 

X Pn+l,n( X ) = Pn+l,n+l( X ) + b 2n+lPn+l,n( X ) + C 2n+lPn,n( X ) + d 2n+lPn,n~l( X ) ■ ( 5 - 15 ) 

By comparing the coefficients of x 2n+1 , x 2n and x 2n ~ l we get expressions for &2n+i? C2n+i 
and d 2n +i respectively in terms of the coefficients a" m and after working out these ex- 
pressions we get the required formulas. ■ 

Observe that the recurrence coefficients have the asymptotic behavior 

b n ~ 3n 2 , c n ~ 3ra 4 , c n ~ ra 6 . (5.16) 
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